Reflectance Fluctuations in an Absorbing Random Waveguide 



T. Sh. Misirpashaev a ' 6 and C. W. J. Beenakker a 
a Instituut-Lorentz, University of Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 
b Landau Institute for Theoretical Physics, 2 Kosygin Street, Moscow 117334, Russia 



Abstract 



We study the statistics of the reflectance (the ratio of reflected and incident 
intensities) of an iV-mode disordered waveguide with weak absorption 7 per 
mean free path. Two distinct regimes are identified. The regime 
shows universal fluctuations. With increasing length L of the waveguide, 
the variance of the reflectance changes from the value 2/15N 2 , characteristic 
for universal conductance fluctuations in disordered wires, to another value 
1/8N 2 , characteristic for chaotic cavities. The weak-localization correction to 
the average reflectance performs a similar crossover from the value 1/3N to 
1/4JV. In the regime 7./V 2 <C 1, the large- L distribution of the reflectance R 
becomes very wide and asymmetric, P(R) oc (1 — R)~ 2 for R <C 1 — jN. 

PACS numbers: 78.20.Ci, 42.25.Bs, 05.40. +j 
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An elegant and fundamental description of universal conductance fluctuations is provided 
by random- matrix theory |]J . Different complex physical systems can be classified into a few 
universality classes, characterized by the dimensionality of the geometry and by the symme- 
tries of the scattering matrix. The so-called circular ensemble, with a uniform distribution 
of the scattering matrix on the unitary group, describes chaotic cavities 0. The variance of 
the conductance in this ensemble (in units of 2e 2 /h) equals 1/8/3, where j3 = 1 (/3 = 2) for 
systems with (without) time-reversal symmetry in the absence of spin-orbit interaction [Q. 
This is the zero-dimensional limit, corresponding to a logarithmic repulsion of the transmis- 
sion and reflection eigenvalues. A disordered wire belongs to a different universality class 
(one-dimensional limit), with a non- logarithmic eigenvalue repulsion and a variance 2/15/3. 
The change from 1/8 to 2/15 is due to a weakening of the repulsion between the smallest 
transmission eigenvalues [|J. 

The optical analogue of universal conductance fluctuations is the appearance of sample- 
to-sample fluctuations in the intensity which is transmitted or reflected by a random medium. 
Universality in this case means that the transmitted or reflected intensity, divided by the 
incident intensity per transverse mode of the medium, fluctuates with a variance which is 
independent of the sample size or the degree of disorder || . It is essential for this universality 
that the incident illumination is diffusive, which means that the incident intensity is equally 
distributed over the N transverse modes of the medium. A new aspect of the optical case 
is the possibility of absorption or amplification of radiation. In Ref. it was shown that 
the distribution of reflection eigenvalues in the limit of an infinitely long waveguide with 
absorption or amplification is the Laguerre ensemble of random-matrix theory. The reflection 
eigenvalues i? n , n = 1,2, . . .N, are the eigenvalues of the matrix product rr^ , where r is 
the N x N reflection matrix of the waveguide. In Ref. the fluctuations in the reflected 
intensity were computed for the case of plane-wave illumination. The purpose of the present 
paper is to consider the case of diffusive illumination, in order to make the connection with 
universal conductance fluctuations. 

We consider the reflection of monochromatic radiation with wavenumber A; by a wave- 
guide of length L and width W. The reflectance R is defined as the ratio of reflected and 
incident intensities. For diffusive illumination it is given by 

i? = iV- 1 ^|r nm | 2 = iV- 1 ^ J R n . (1) 

n,m n 

The L-dependence of the distribution of reflection eigenvalues is described by the Fokker- 
Planck equation || 

r\ r\ TV f) 

BP 1 i 
HT+PPT, 7 r + 7(/?iV+2-/3)P , (2) 



X 



where we use the parametrization R n = A n /(1 + A„), A n G (0, oo). The parameter 7 > is 
the ratio of the mean free path I and the absorption length l a . Eq. (0) is valid if kl ^> 1, 
kl a 3> 1, and W L. Optical systems normally possess time-reversal symmetry (/3 = 1), 
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but in view of a recently observed magneto-optical effect [0], we also include the case of 
broken time- reversal symmetry ((3 = 2). 

Relevant length scales are defined in terms of the transmitted intensity |||9| . The trans- 
mitted intensity decays linearly with L for L <C £ and exponentially for L ^> £, where the 
decay length £ = (l/£ a + l/^) -1 contains a contribution £ a = /(27 + 7 2 ) -1 / 2 from absorption 
and a contribution = A((5N-\-2—(3) from localization by disorder. We will study the two 
regimes (1) £ a <C & (or 7iV 2 ^> 1) and (2) ^ a 3> (or 7 A/ -2 C 1). In both regimes we 
assume N 1. 

(1) The regime 7iV 2 ^> 1 admits a perturbative treatment for all L by the method of 
moments of Mello and Stone []!]]. We define the moments M q = A r_1 (— l) q £^(1 + Aj) 



so that i? = 1 + Mi. The Fokker-Planck equation (|j) enables us to express derivatives 
d(M q ^ . . . M™ fe fc )/ <9L in terms of the moments themselves. Expanding each moment in powers 
of 1/N, we get a closed system of first order differential equations for the coefficients. The 
method has been explained in detail |TD[, therefore we just present the result. We need the 



following terms in the 1/iV-expansion of the moments, 

(M?) = F qfl + N- l F q + + N~ 2 F q , 2 + 0(N- 3 ), (3a) 

(M?M 2 ) = G qfi + N- l G q + + 0(N~ 2 ) 1 (3b) 

(M%M 3 ) = H q + 0(N- 1 ). (3c) 

The g-dependence of the coefficients is given by 

F q , = f\ F q>l = 5 m nP~\ F 9 , 2 = q^p- 1 + q(q - l)^/ 9 " 2 , (4a) 

G q ,o = gf q , G q>1 = Spyiaf* + q M f q - X ), (4b) 

H q = hp. (4c) 

Here /, g, h, fi, # lj2 , a, rjx^ are functions of L which obey the following system of first order 
differential equations and initial conditions: 



Idf/dL = 


f- 


27/ -2 7 , /(0) = -l, 






(5a) 


Idg/dL = 


4(/- 


-7k + 2(/ 2 -2 7 /), 5 (0) = 1, 






(5b) 


Idh/dL = 


6(/- 


- 7 )h + 3g(g + 2f-2 7 ), h(0) = -l, 






(5c) 


Idfi/dL = 


2(/- 


-j)fi + g-f 2 , M0) = 0, 






(5d) 


IdOp/dL = 


= 4(/ 


- 1 )9 p + (2/P)(g + h)+6p lf i(9-f 2 ), 


W) 


= 0, 


(5e) 


Ida/dL = 


4(/- 


- 7)0- + 4/iQ? + / - 7) + % - f) + 4(/i - 


/<?), 


a(0)=0, 


(5f) 


Idrjp/dL - 


= 2(f 


- 7)7^ + 2^ + 5^ (or - 2nf - g + / 2 ), 


^(0) 


= 0. 


(5g) 



The equations can be easily solved one after another. The first two of them determine the 
mean and variance of the reflectance for plane- wave illumination, studied in Ref. ||. The 
other equations determine the mean and the root-mean-squared fluctuations for diffusive 
illumination to order N" 1 . We decompose (R) = Ro + SR, where Ro = O(N ) and 5R = 
0(N~ 1 ). In terms of the functions in Eq. fl5|) we have 

R = l + f, 5R = Sfnfi/N, (6a) 
v&rR = (26 p - 5 pi ^)/N 2 = A6 2 /{3N 2 . (6b) 



3 



0.135 

cr 0.130 

ni 0.125 

: 0.120 

0.115 

0.30 
So 0.25 

z 

0.20 
0.15 



Z 1 1 I I I I I l| 1 1 I I I I I l| 1 1 I I I I M 




J I I I I 



1 1 I I I I I l| 1 1 I I I I I l| 1 1 I I I I M 



J i/ I I I I 



10 100 1000 
L/Z 



FIG. 1. Length dependence of N5R and N 2 vaiR for (3 = 1 and 7 = 10 -4 , according to 
Eqs. (H), (|6j), in the regime 7./V 2 S> 1. The variance vai R of the reflectance crosses over from a 
plateau at 2/lhN 2 (one-dimensional limit) to a plateau at 1/8N 2 (zero-dimensional limit). The 
crossover is non-monotonic and occurs when the length of the waveguide becomes comparable to 
the decay length £ = 70 1 of the transmittance. The weak-localization correction 5R shows a similar 
crossover from 1/3ZV to 1/4ZV, but the plateaus are less well-defined for this value of 7. 

The weak-localization correction SR vanishes for (3 = 2, while vari? oc 1/(3, just as without 
absorption. 

To find (R) and var R we need to solve the first five equations (g). The analytical solution 
is cumbersome to use, but it is easy to integrate the system numerically at a given 7. Results 
for 7 = 10~ 4 are shown in Fig. |l[ The large-L limit can be found directly by replacing the 
derivatives at the left-hand side of Eq. (Q) by zero. In this way, we obtain the following 
asymptotical values: 



Ro ■ 
5R 

var R 



1 +7 - ^27 + 72 



SpiN- 



1 (K 2 + 7)^ + 7 " (1 + 7)V / 7(2 + 7)- 1 



2(3N 2 (2 + -ff 



(7a) 
(7b) 

(7c) 



The large-L limit (|7|) can also be obtained from the Laguerre ensemble for the reflection 
eigenvalues 0, 
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Poo({A n }) cc exp [p ^ I A, - \\ - 1 ((3N+2-(3) ]T A,] , (8) 



which is the asymptotic solution of the Fokker-Planck equation (g). The limiting values (0) 
are reached when L>(, hence when the transmittance through the waveguide has become 
exponentially small. For L <C £ the effect of absorption can be neglected. Over a range of 
lengths L such that I C L <C the mean and variance are given by QTOJ i?o = 1 — l/L, 



SR = 6pi/3N, vaiR = 2/15(3 N 2 . The value 1/8/? iV 2 for the variance when 7 < 1 and 
L ^> £ follows directly from the logarithmic repulsion of the A n 's in the Laguerre ensemble 
Ijnj . The difference with the value 2/15/5 N 2 for I <C L <C £ arises because the repulsion is 
non- logarithmic in the absence of absorption @. 

^ We now turn to the second regime, 'jN 2 <C 1. This regime cannot be treated 
perturbatively for L > £, because of the onset of localization by disorder. The limiting 



large- L values of (R) and vari? can be computed from Eq. ([8]) using formulas |12| for the 



density and correlation function of the A n 's in the Laguerre ensemble. The result is 

(R) = l- PlN [ln(l/ 7 iV 2 ) + Oil)] , (9a) 
var R = (3j- {(3-fN) 2 [ln(l/ 7 iV 2 ) + 0(1)] 2 . (9b) 

A new crossover length L c = £ ln(l/7iV 2 ) appears at which the mean and variance of R attain 
their asymptotic values @. This is the length at which the exponential decay (1 — R) ~ 
N' 1 exp(— L/4£) comes to a halt. 



From Eq. (§) we see that the reflectance has a wide distribution in the regime 7 iV 2 <C 1: 
The root-mean-squared fluctuations of 1— R are greater than the mean by a factor ^iV 2 )" 1 / 2 . 
(In the opposite regime 7iV 2 ^> 1, in contrast, the distribution of the reflectance is a narrow 
Gaussian, see inset in Fig. |2|.) To determine the tail of the distribution, it is sufficient 
to consider only the contribution from the smallest eigenvalue Ai, which gives the main 
contribution to 1 — R when 7iV 2 — > 0. The smallest eigenvalue in the Laguerre ensemble 
has the exponential distribution P(Ai) = ft^N 2 exp(— fl^N 2 ^) [pE|, hence 



P(R) = p-yN(l - R)~\ -fN < 1 - R < l/N. (10) 

We have calculated P(R) by generating a large number of random matrices in the Laguerre 
ensemble (see Fig. 0). The distribution reaches its maximum at 1 — R ~ 7^ and then drops 
to zero for smaller values of 1 — R. The tail for large values of 1 — R is well described by 
Eq. ( [Top (solid curve in Fig. 0). 

To conclude, we have studied the statistics of the reflectance in an absorbing random 
waveguide under diffusive illumination. When the decay of the transmittance is dominated 
by absorption ( 7 iV 2 ^> 1), the fluctuations are shown to possess the same features as univer- 
sal conductance fluctuations, including independence on the disorder and 1/(3 dependence 
on the symmetry index. A crossover from zero- dimensional to one-dimensional limit was 
found in Refs. |I5| , |T6|| for a chain of chaotic cavities. (A long chain of cavities behaves as 



a diffusive wire.) We have found an opposite crossover, from the one-dimensional to the 
zero-dimensional limit, as the length of the waveguide is increased beyond the decay length 
£. Another regime, when the decay of the transmittance is dominated by localization due 
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FIG. 2. Large-L distribution of the reflectance for (3 = 1. Data points are obtained by generat- 
ing 5-10 4 random matrices in the Laguerre ensemble The main plot is for the regime 7./V 2 <C 1 
(N = 20, 7 = 10 -4 ), the solid curve being the asymptotic tail (|l0|). The inset is for the regime 
■yN 2 ^> 1 [N = 20, 7 = 0.1), the solid curve being a Gaussian with mean and variance given by 
Eq. (0). 

to disorder {^N 2 <C 1), is principally new and characterized by a wide and asymmetric 
distribution of the reflectance. The asymptotic regime establishes at a new characteristic 
scale L c > £. 
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